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A  METHOD  FOR  STUDYING  THE  HUNTING  OSCILLATIONS  OF  AN  AIRPLANE  WITH  A  SIMPLE 

TYPE  OF  AUTOMATIC  CONTROL 

By  Robert  T.  Jones 


SUMMARY 

/I  method  is  presented  for  predicting  the  amplitude  and 
frequency,  under  certain  simplifying  conditions,  of  the  hunting 
oscillations  of  an  automatically  controlled  aircraft  with  lag  in 
the  control  system  or  in  the  response  of  the  aircraft  to  the  controls. 
If  the  steering  device  is  actuated  by  a  simple  right-left  type  of 
signal,  the  series  of  alternating  fixed-amplitude  signals  occurring 
during  the  hunting  may  ordinarily  be  represented  by  a  “square 
wave."  Formulas  are  given  expressing  the  response  to  such  a 
variation  of  signal  in  terms  of  the  response  to  a  unit  signal. 
A  more  complex  type  of  hunting,  which  may  involve  cyclic 
repetition  of  signals  of  varying  duration,  has  not  been  treated 
and  requires  further  analysis.  Several  examples  of  application 
of  the  method  are  included  and  the  results  discussed. 

INTRODUCTION 

When  an  airplane  or  otlior  aircraft  is  directed  by  a  simple 
right-left  signal  from  an  automatic  steering  device,  the 
result  is  usually  a  maintained  hunting  oscillation  about  the 
desired  path.  The  amplitude  of  this  oscillation  is  influenced 
by  the  amount  of  backlash  or  “dead  spot”  in  the  control 
system  and  by  the  damping  of  the  motion  of  the  airplane. 
In  the  following  analysis  the  amplitude  and  frequency  of 
these  oscillations  is  investigated  in  terms  of  the  response 
characteristics  of  the  airplane. 

ANALYSIS 

The  analysis  is  based  on  consideration  of  the  response  of 
the  airplane  (in  terms  of  angle  of  yaw  or  pitch)  to  a  continued 
(unit)  signal  (fig.  1).  This  response  may  be  calculated  by 
the  ordinary  theory  of  dynamical  stability  and  is  conveni¬ 
ently  represented  in  operational  form  (references  1  and  2)  as 
follows: 

Ri(f)=Ri  (.0)1(0  (1) 

The  unit  response  ordinarily  occurs  in  the  form 

Si  \D)  p(D) 

from  which  is  obtained 

Rft)= <7(0  -K<V“+<V*‘+  •  •  •)  (2) 

where  C(t)  is  the  steady-state  motion,  and  Ch  are  the 
constant  coefficients  of  the  Heaviside  expansion,  and  Xj,  X2, 
and  so  forth,  are  the  nonzero  roots  of  the  characteristic 


equation  defining  the  natural  periods  of  oscillation  and  the 
damping  of  the  aircraft  without  signal.  The  function  f(D) 
and  the  particular  solution  (7(f)  depend  on  the  time  variation 
of  control  displacement  produced  by  a  signal  and  on  the 
stability  characteristics  of  the  airplane  in  the  degrees  of 
freedom  in  which  the  control  operates.  (See  reference  3.) 
In  the  case  of  a  continued  signal,  the  usual  form  of  the  func¬ 
tion  (7(f)  is 

(7(f) =<7-1-1- Oof 

where  (70  is  the  steady  rate  of  turn  called  for  by  the  signal 
and  O-i  is  a  constant.  (See  fig.  1.) 


Figure  1. — Typical  response  to  continued  signal. 

During  a  hunting  oscillation,  the  automatic  steering  device 
reverses  the  signal  periodically  as  the  airplane  swings  through 
the  desired  heading.  A  typical  hunting  oscillation  is  shown 
in  figure  2.  Here  it  is  assumed  that  the  reversal  of  signal 
is  delayed  either  because  of  a  “dead  spot”  in  the  steering 
device  or  because  of  backlash  in  the  control  mechanism  or  a 
combination  of  the  two.  As  indicated,  the  oscillation  will 
have  a  fundamental  period  2r/co  (where  a  is  the  angular 
frequency  of  the  hunting  oscillation)  but  may  also  involve 
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components  of  higher  frequency,  depending  on  the  natural 
modes  of  oscillation  of  the  airplane.  Ordinarily,  the  shorter- 
period  components  do  not  have  sufficient  amplitude  to  cause 
a  reversal  of  the  signal  during  a  half  cycle.  In  these  cases 
the  variation  of  signal  with  time  will  be  represented  by  a 
simple  “square  wave,”  which  may  be  expressed  as  a  function 
of  time  by 

-  si11  nut  (n= 1,  3,  5,  .  .  .)  .  (3) 

■tr  „  n 


dead  spot.  With  the  frequency  determined,  it  is  possible 
also  to  find  the  amplitude  of  the  oscillation  and  the  maximum 
deviation  of  the  airplane  from  its  path. 

In  the  simplest  cases  the  required  information  may  be 
obtained  directly  from  equation  (6).  In  the  case  of  more 
complex  motions,  further  analysis  will  be  required  as  follows: 

As  a  first  step,  separate  Ri(iruo)  into  its  real  and  imaginary 
parts 

Ri  (incS)  =A(ruo)  -j-iB  (nai) 


or,  more  convenieutly,  by  the  imaginary  part  of  the  cor¬ 
responding  exponential  series;  that  is, 


The  functions  A  and  B  may  be  plotted  against  nu  as  in 
figure  3.  These  functions  will  show  peaks  near  values  of  nu 


4  1 

I.  P.  -2-  «*“* 

■KnU 


(4) 


where  is  taken  to  represent  a  time  at  which  the  signal 
becomes  positive. 


The  response  to  the  alternating  signal  is  obtained  by 
substituting  expression  (4)  for  the  unit  function  1  (£)  in 
equation  (1).  Thus, 

Rit^I.P.RdD)-^#'*  (5) 

Tf  n  n 

If  the  airplane  is  inherently  stable,  so  that  transient  effects 
following  the  start  of  an  oscillation  disappear  with  time,  the 
remaining  steady  oscillation  will  be  represented  by 


Figure  3.— Curves  showing  In-phasc  and  oat*of-phose  components  of  response  to 
periodic  signal. 


corresponding  to  the  resonant  frequencies  of  the  airplane. 
Then,  for  any  particular  hunting  frequency  o>, 

R(t)=—  y1.—  [A{ru o)  sin  nwt-\-B(nw)  cos  nui]  (7) 

7T  n  71 

At  the  time  of  reversal  of  the  signal  sin  no>t=0  and 
cos  jwt=  ±  1,  the  sign  depending  on  whether  the  signal  is 
becoming  positive  or  negative.  The  amplitude  of  the  response 
at  this  instant  is  therefore 


R(t)=I.  P.  -  y~  Rl(j,nw)ei*'“  (6) 

TT  n  n 

Equation  (6)  gives  the  forced  oscillation  of  the  airplane  in 
response  to  an  alternating  signal  in  the  form  of  a  square  wave 
of  any  frequency  a. 

By  investigating  the  form  of  these  forced  oscillations  at 
various  frequencies  it  will  be  possible  to  ascertain  whether 
such  oscillations,  under  the  conditions  of  automatic  control, 
will  give  rise  to  the  assumed  alternating  signals  of  equal 
duration,  and  thus  to  establish  certain  ranges  of  u  over  which 
hunting  of  this  type  can  occur.  It  will  also  be  possible  to 
establish,  in  these  ranges,  a  correspondence  between  the  fre¬ 
quency  of  the  hunting  oscillation  and  the  magnitude  of  the 


±1\b(.u)+\b(Zo>)+\b{50>)+.  .  .] 

This  amplitude  will  also  be  the  amplitude  of  the  dead  spot. 
(See  fig.  2.)  A  plot  of 

(n—1,  3,  5,  .  .  .) 

7 1  n 

can  readily  be  obtained  from  the  curve  of  B  in  figure  3  and 
will  show  the  periods  of  the  hunting  oscillation  corresponding 
to  various  widths  of  dead  spot. 

The  slope  of  the  response  curve  at  this  same  instant  is 

-S/s=^^)b=—  [A(w)+A(3co)+A(5w)-f-  .  .  .] 
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If  the  response  to  a  positive  signal  is  negative  (as  in  fig.  1), 
in  order  that  the  motion  represent  a  possible  hunting  oscilla¬ 
tion  (that  is,  be  consistent  with  the  assumed  variation  of 
signal),  it  is  necessary  that 

( I )  RB  ^  0 
for  a  positive  dead  spot,  and  that 

(II)  R'b>  0 


EXAMPLES 

In  order  to  demonstrate  and  check  the  procedure  described, 
assume  a  simple  response  characteristic  in  which  the  air¬ 
plane  immediately  starts  turning  at  a  constant  rate,  as 
directed  by  the  signal.  With  this  response 


indicating  that  the  airplane  crosses  the  dead  spot  in  the 
proper  direction.  A  further  condition  is  that  no  more  than 
one  complete  crossing  of  the  dead  spot  occurs  within  one- 
half  cycle;  that  is, 

(III)  R{t)>-R0 

(See  fig.  2.)  The  value  of  R(t)  in  the  middle  of  a  half  cycle 
is  relatively  simple  to  obtain  '( 

Ja(u)  ~A(3o>)  +|A(5co)  -  .  .  .] 


and  may  be  used  as  a  criterion,  though  RA  is  not  necessarily 
the  maximum  or  minimum  value  of  R{t)  (see  fig.  4)  and 
condition  III  may  not  be  satisfied  even  though  RAy>—RB. 


Figure  4. — Plot  of  Ra  and  Rb  against  frequency,  showing  approximate  regions  in  which 
hunting  oscillations  aro  possible  and  width  of  dead  spot  in  those  regions. 

,  It  should  be  noted  that,  in  the  regions  excluded  by  the 
foregoing  conditions,  a  more  complex  type  of  hunting 
oscillation  involving  a  sequence  of  signals  of  different  dura¬ 
tions  may  occur.  In  these  regions,  the  curves  of  RA  and 
Rb  derived  for  the  square-wave  signal  no  longer  apply  to 
the  condition  of  automatic  control.  These  oscillations  re¬ 
quire  analysis  beyond  that  presented  in  this  report. 


*1(7icof)=.— °i 
'  TUO 

and,  from  equation  (7), 

R(t) =-S  ~T  cos  Aw* 
irT  n2a 

which  is  the  Fourier  series  for  a  “saw-tooth”  wave  90°  out 
of  phase  with  the  signal.  (See  fig.  5.)  In  this  case  the 
response  occurs  without  lag  and  the  amplitude  of  the  hunt¬ 
ing  is  exactly  equal  to  the  dead  spot.  The  frequency  t» 
is  tCo/2  divided  by  the  width  of  the  dead  spot. 


A  simple  example  nearer  the  practical  case  is  one  in- which 
the  signal  causes  a  force  F  to  act  on  a  mass  m.  In  this  case 
the  response  to  a  unit  signal  is 


and  the  hunting  oscillation  is  seen  to  be 

R  (t)  =~  sin  nwt 

7T  m „  n?or 


The  expression  is  recognized  as  the  Fourier  series  for  a 
succession  of  parabolic  segments  (fig.  6).  It  should  bo 
noted  that  there  is  no  component  out  of  phase  with  the 
signal,  with  the  result  that  RB  is  zero  for  all  values  of  to_ 
Hence  the  calculation  shows  no  possibility  of  hunting  withj 
a  finite  dead  spot.  In  fact,  it  can  be  seen  from  energy 
considerations  that  if  a  dead  spot  existed  the  oscillation 
would  be  divergent. 
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Figure  6. — Hunting  oscillation  of  mass  nctod  on  by  force. 


Interesting  applications  of  the  method  are  furnished  by 
cases  in  which  the  response  to  a  signal  shows  a  lag,  possibly 
■due  to  backlash  in  the  control  mechanism,  in  addition  to  a 
dead  spot.  A  simple  example  of  this  kind  is  illustrated  in 
figure  7.  Here  the  response  is  similar  to  that  in  the  first 
example  (fig.  5)  except  for  the  time  lag  t.  Use  is  made  of 
the  well-known  lag  operator  e ~rD.  Thus, 

e-rom=j{t_r) 

Applying  this  operator  to  the  response  in  figure  5  gives 


Rl(D)  =  -e 


R.  (inw)  =—  (sin  ncor+i  cos  nar) 

'  7 10} 


and,  finally  (equation  (7)), 


/ 

/ 


1  / 

Figure  7.— Example  In  which  response  shows  lag. 

reference  to  figure  7  that  the  half  period  of  the  oscillation 
in  this  case  (no  dead  spot)  is 


0) 


Langley  Memorial  Aeronautical  Laboratory, 
National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  May  5,  1944- 


4  <7  . 

R(t)=—  yy_S-  (gin  7 noT  sin  «co<+ cos  i uor  cos  na it) 
v '  v  *  n 


=—  cos  7Uo(<  — t) 

■K  „nra 


With  the  lagging  response,  the  hunting  oscillation  is  not 
confined  to  the  amplitude  of  the  dead  spot  and,  in  fact, 
hunting  will  occur  with  no  dead  spot.  It  is  easily  seen  by 
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